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Введение. Исследуется влияние двух способов организации скалярного произведения на скорость сходимости 
решения в энергетическом методе граничных состояний. Основу метода исследования составляют пространства 
внутренних и граничных состояний, которые сопряжены изоморфизмом. Оба пространства ортонормируются, 
используя то или иное скалярное произведение. Искомое состояние раскладывается в ряд Фурье по элементам 
ортонормированного базиса, определяются коэффициенты этой линейной комбинации. Различие двух способов 
заключается в назначении скалярных произведений и вычислении коэффициентов Фурье. 

Материалы и методы. Применительно к методу граничных состояний предложена новая теория организации 
скалярного произведения в пространствах внутренних и граничных состояний. Построены вычислительные 
алгоритмы ее практической реализации. В традиционном (первом) способе в качестве ортогонализатора в про- 
странстве внутренних состояний используется внутренняя энергия упругого деформирования. Здесь коэффици- 
енты Фурье представляют собой работу заданных сил на базисных векторах перемещения точек границы. В 
исследуемом (втором) способе скалярные произведения представляют собой интегралы от перекрестных про- 
изведений базисных векторов сил на границе. Соответственно коэффициенты Фурье вычисляются как интегра- 
лы произведения заданных сил на границе тела на базисные векторы сил. 

Результаты исследования. Проведено численное исследование первой основной осесимметричной задачи тео- 
рии упругости для трансверсально-изотропного цилиндра при отсутствии и при наличии массовых сил. При 
отсутствии массовых сил анализ упругих полей, полученных при одинаковом числе используемых базисных 
элементов, показал, что второй способ имеет наибольшую точность результатов. При наличии массовых сил 
второй способ не показал эффективности в плане единственности решения, однако он вполне пригоден для по- 
строения множества упругих полей, используемых в решении более сложных задач. 

Обсуждение и заключения. Полученные результаты могут быть использованы в решении краевых задач меха- 
ники не только анизотропного, но и изотропного тела. При решении более сложных задач, таких как контакт- 
ные и смешенные, вопрос о скорости сходимости требует отдельного исследования. 
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стояний, первая основная задача, массовые силы. 


Образец для иитирования: Иванычев, Д. А. Два способа организации скалярного произведения в методе гра- 
ничных состояний / Д. А. Иванычев // Вестник Донского государственного технического университета. — 2020. 
— Т. 20, № 1. — С.15—24. 6 рз://4о1.оге/10.23947/1992-5980-2020-20-1-15-24 


Финансирование: Исследование выполнено при финансовой поддержке РФФИ и Липецкой области в рамках научного 
проекта № 19-41-480003 "ра". 


© Иванычев Д. А., 2020 





Тууо ууау$ оЁ огоап171то зсааг ргодисЕ ш Пе Боипдагу $6ж{е те од 


О. А. ГуапусВеу 
Глрезк Зе Тесбиса| Олтуегзиу (ГалребК, Визз1ап Еедеганоп) 


[итодисноп. Те шЙоепсе о{ {\о \ауз оЁ огоап1ттае зса[аг ргодис{ оп Ше сопуегоепсе ге оЁ Ше зошНоп ш Ще еп- 
егоу тео о? Боппдагу з1а{ез 15 сопз14еге4. ТВе тео 1$ Базе оп Ше 5расез оЁ пщегпа| ап Бочпдагу зез \мсВ 
аге сопласайеа @гочов 15о1тогр1$т. Во 5расез аге огопогта|П7е4 изо опе зса1аг ргодис{ ог апо!ег. ТБе дезшед 


Механика 


15 


БИр://уезыиК.Аоп$а.га 


Вестник Донского государственного технического университета. 2020. Т. 20, № 1. С. 15-24. 155М 1992-5950 е155М 1992-6006 
Тези! о{ оп ие Тесйтса От»егхйу. 2020. Иа. 20, по. 1, рр. 15—24. 155М 1992-5980 е155М 1992-6006 





убае 15 ехрапае ш Фе Еочпег земез ассог4ше фю Фе еетеп5 о? Фе обропогта|7е4 Ъа$1$; ап бе соеЁЯслет$ оЁ 
{61$ Ппеаг сотЫпайоп аге д&егиитеа. ТБе мо шефоа5 ег ш Фе азз1етитепЕ оЁ зсайаг ргод4ис{ ап Фе сайсиайоп 
ог фе Еоипег соеЁЙслепт 5. 

Маептай апа Мефос$. Шш те|аНоп о фе тефо4 о? Ъоипдагу ${аез, а пе\ми Феогу оЁ огоап1ттае а зсайаг ргодиасе ш фе 
зрасез оЁ ищегпа| ап Боипдагу зе$ 15 ргорозе4. СотрщаНопа| а1оогИ из аге сопзгаце Юг {$ ргасйса| паре- 
тетаноп. ш Фе гад1опа| (Вг$ё) арргоасй, ®е пуегпа[ епегоу оЁ еазис аеГогтайоп 1$ изе4 аз ап оппогопайтег ш 
Фе зрасе ог пщегпа| аз. Неге, фе Еочиег соеЁЯ слеп аге Фе ууогК оЁ о1уеп Ююгсез оп Фе Баз1$ уесфогз ог 41зр1асе- 
еп о? Ще Боцпдагу рош. ш Фе зад1е4 (зесопа) арргоасВ, зсайаг ргодис 5 аге ищеота1$ ог Ше сгоз$ ргодис{$ ог е 
Ба$1$ Югсе уесогз ай фе Боцпдагу. Ассогаше]у, Фе Еоцпег соеЁИслет5 аге са1ся]=ме4 аз пеога|$ оЁ Ше рго4дис{ оЁ 
Фе отуеп Югсез а фе Боду Боцидагу Бу Фе Баз1с Гогсе уесюгс. 

Кези5. А питепса1 заду оЁ Фе Вгз{ ргитла[ ах1зутлтейлс рго ет оЁ Фе е]азисйу Шеогу Юг а Напзуегзе]у 1зогор1с 
суппдег ш Фе абзепсе апа ргезепсе о{ тпаз$ Ююгсез 15 сопдицеа. ш Фе абзепсе оЁ та$$ Гогсез, ап апа[у$1$ о? Ве еаз- 
ис Не4$ оМашеа Гог Фе зате питьфег оЁизеа Баз1с @етет Ваз зво\уй фа{ Ще зесоп4 тео4 Ваз Ве отезжез{ асси- 
гасу оЁ Фе гезиИз. Оп4ег зо]ушя Фе ргоет у Ше ргезепсе ог таз юогсез, фе зесоп4 тефо4 414 по{ зВо\и еЁЙ- 
слепсу ш {егил$ оЁ Фе итацепе$$ оГ Фе зоаНоп; Во\еуег, И 15 даце зайца е Юг сопзгасйипе а тиаде о? @аз$ис 
Не!4$ изе4 {0 зо]уе тоге сотр] ех ргое!т. 

Оксиуятоп апа Сопсияоп5. ТВе гези{$ оМаше4 сап Бе изе4 ю зо]уе Боипдагу-уаше рго]етаз о? тесВап1с$ оЁ по 
оту ап ап1зогор1с Боду, Би а[50 ап 15 огор1с опе. УМВеп 50!уше тоге сошр!ех ргоептл$, засВ аз сощас{ апа пихе4 
опез, Ве 155ие оР Пе гже оЁ сопуегоепсе гедиигез а зерагае заду. 


Кеуногау: Боипдагу $айе тефо4, зсаЙаг ргодис®, ищегпа] епегоу, %{ае зрасез, {Ве Иг$Е таш фазК, таз$ Гогсез. 


Еог сйапоп: О. А. [уапусВеу. Туо \ауз оЁ огоап1тлае зса|аг ргодисе ш фе Боцпдагу зе те од. УезникК о ОТО, 
2020, уо1. 20, по. 1, рр. 15-24. В рз://401.0те/10.23947/1992-5980-2020-20-1-15-24 


Еипат> трогтаноп: Тфе заду \аз сагле4 ошё уп Фе Нпапса| заррог о ВЕВЁК апа Фе ГлреёзКк Веслоп аз рам оЁ 
Бе гезеагсВ рго]ес+ М. 19-41-480003 "р_а". 


Введение. Краевые задачи теории упругости в механике достаточно изучены, поэтому в последние го- 
ды проводятся исследования частного направления. Например, рассмотрена осесимметричная первая основная 
задача для полуполосы [1], решение которой строится в виде разложений по системам функций Фадля- 
Папковича и имеет явный вид. Предлагается общий метод решения первой основной задачи теории упругости 
для прямолинейно-анизотропного тела в случае плоской деформации [2]. Используются замкнутые системы 
краевых задач, схожих с задачами Гильберта, что позволяет добиться большей общности метода. Методом ко- 
нечных элементов на основе вариационного принципа Кастильяно решены плоские изотропные задачи [3]. Это 
позволило получить поля напряжений на сетках достаточно низкой размерности, в том числе для несжимаемых 
материалов. Решаются контактные задачи о внедрении эллиптических штампов в трансверсально-изотропное 
упругое полупространство [4]. 

Метод граничных состояний при решении краевых задач для анизотропных тел также нашел свое при- 
менение. Например, рассмотрено упругое равновесие трансверсально-изотропного цилиндра под действием 
осесимметричных поверхностных сил [5]. Исследуются задачи кручения протяженных цилиндров из материала 
с анизотропией общего вида [6]. Показана математическая модель получения явных параметрических решений 
для изотропных и анизотропных тел [7, 8], где константы среды в качестве параметров включены в упругие 
поля. Разработана методика решения задач теории упругости средствами компьютерной алгебры [9]. В реше- 
нии задач напряженно-деформированной неограниченной упругой среды, содержащей сферические полости 
или включения, при разных условиях применен метод граничных состояний [10]. 

В настоящей работе исследуются два подхода к назначению скалярного произведения в «теле» метода 
граничных состояний. При этом проводится тестирование каждого состояния на примере решения первой ос- 
новной задачи теории упругости. В каждой задаче удерживается одинаковое число используемых элементов и 
оценивается уровень погрешности. 

Материалы и методы. Метод граничных состояний (МГС) [11] является энергетическим, в нем при- 
меняется фундаментальная теория рядов для решения основных задач механики. В качестве опорных исполь- 
зуются понятия внутренних и граничных состояний. Внутреннее состояние & обусловлено набором вектора 


перемещений и, тензора деформаций = и тензора напряжений Т: 
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6={м,в,Т}. (1) 
Граничное состояние обусловлено набором вектора перемещений точек границы и’ и усилий р на границе те- 
ла: 


__ у 
у=\и ‚РН 
Совокупность таких состояний образует базисы пространств внутренних = = ео и граничных 
Г = Ре. состояний. Далее проводится ортогонализация базисов состояний, где в качестве орто- 
гонализатора в базисе граничных состояний используется выражение: 
(5,,5,)=Г а ТАУ, 
7 
в базисе граничных состояний — выражение: 
__ у 
(уру,) = р,а5. (2) 
5 


Каждому элементу 6, е= соответствует единственный элемент 7, Е /`’, причем это соответствие взаимно- 
однозначное: &, <> у, . Это позволяет отыскание внутреннего состояния свести к построению изоморфного ему 


граничного состояния. В случае первой основной задачи искомые внутреннее и граничное состояния представ- 
ляет собой ряд Фурье: 


[е.9) [®.9) 
=> с, => СТ, , (3) 
К=1 К=1 
здесь с, — коэффициенты Фурье: 
с; = | ри, 45, (4) 
К) 
где р — вектор заданных поверхностных сил; и, — вектор перемещения в К-ом базисном элементе базиса 


внутренних состояний. При этом базисные наборы формируются на основании общего или фундаментального 
решения задачи. 

Первые общие решения уравнения Ламе линейной теории упругости были построены еще в 30-е годы 
прошлого столетия. Уравнение Ламе — это уравнение движения (в рассматриваемом случае — равновесия) 
Эйлера: 

УТ+ Г =0, 
где Т — тензор напряжений; У — оператор Гамильтона, действующий как дивергенция; } — массовые силы. 


В уравнении Ламе тензор напряжений Т в соответствии с законом Гука представлен через тензор де- 
формаций =. В свою очередь тензор деформаций в соответствии с соотношением Коши представлен через век- 
тор перемещений и. В общих решениях уравнения Ламе вектор перемещения определяется: в теории изотроп- 
ной упругости через гармонический вектор В и гармонический скаляр, в теории анизотропной упругости — 


через функцию напряжений Е. 


Вектор В (функция Е) может быть представлен в виде ряда по базисным векторам В, = В, (в) — 


функциям координат а’. В результате каждому гармоническому базисному вектору В, (функции Е.) будут 
поставлены в соответствие следующие базисные элементы: 

» вектор перемещения и, ; 

»® тензор деформаций &, ; 

® тензор напряжений Г, ; 

»® вектор массовых сил /, (из уравнений равновесия); 

» вектор поверхностных сил (из фундаментального соотношения Коши): 

р, =П1,, 

где и — внешняя единичная нормаль к поверхности тела. 


По перечисленным базисным элементам соответствующие векторы или тензоры разлагаются в ряды 
Фурье с одинаковыми коэффициентами с,, которые определяются из условий ортогональности базисных 


функций. Например, для первой основной задачи при отсутствии массовых сил, когда на всей поверхности тела 
5 заданы внешние силы р и осуществлена ортогонализация базисных векторов р, (|р.р ‚45 =6, — дельта Кро- 
К) 


некера), коэффициенты с, определяются из выражения: 
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с =|рр,95. (5) 
5 
Это выражение вытекает из представления: 
р=>.сРь. 
К=1 


Таким образом, исследуется способ формирования решения с помощью выражения для скалярных 
произведений в базисе граничных состояний: 


(у›7,) =] р.р5 (6) 


и выражения для коэффициентов Фурье (5). 
В случае второй основной задачи имеют место зависимости: 


(у›у,) = [ги 45; 
К) 
с =[ии, 5, 
К) 


где и — заданный вектор перемещения точек границы тела; и, — вектор перемещения в А-ом базисном эле- 


менте базиса граничных состояний. 

Результаты исследования. Параметры скорости сходимости рядов и точности результата рассмотрим 
на примере решения задачи об упругом равновесии трансверсально-изотропного цилиндра из темно-серого 
алевролита [12] в безразмерном виде (рис. 1). 

Граничные условия: 


р=О0, 5 2=—-2, 0<и<1; 
р=О0, 5, | 2=2,0<и<1; 


р, =4-2,р, =0,5,| г=,-2<252. 


ее — 





Рис. 1. Граничные условия для транстропного цилиндра 


Подробно методика решения первой основной задачи традиционным способом с помощью скалярного 
произведения (2) приведена в работе [6]. Базис внутренних состояний (1) строится следующим образом: 

® используя общее решение задачи о плоской деформации [13], конструируются базисные наборы 
плоских вспомогательных состояний; 

» по формулам перехода определяются базисные наборы пространственных осесимметричных состо- 
яНий; 

е проводится ортонормирование базиса внутренних состояний по матричному алгоритму Грама- 
Шмидта с использованием скалярного произведения (2); 

® из ортонормированного базиса внутренних состояний редуцируется ортонормированный базис гра- 
ничных состояний; 

е вычисляются коэффициенты Фурье (4) и строятся ряды (3) в развернутом виде (индекс К помещен 
наверх): 


[©.0) [©.0) [©.0) 
Е. —_ к. —_ к. — р 
Си; › Р; 22. СьР Бо; = 2.19 › 5; = 2. 48, 


в 


и. = 


1 


К 


| 


1 
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Опустим информацию о полях характеристик напряженно-деформированного состояния, полученных 
при одном и другом способе назначения скалярного произведения, и приведем лишь основные результаты. При 
этом назовем традиционный подход, используемый в [11], как первый способ решения задачи, а подход, ис- 
пользующий скалярное произведение (6) и коэффициенты Фурье (5), как второй способ. Точность решения при 
удержании одинакового числа базисных элементов у второго способа выше. На рис. 2 для каждого способа 
приведено сопоставление полученных граничных условий (ГУ) с заданными при использовании 8-ми значений 
коэффициента Фурье. Усилия изображены в масштабе, например, истинное значение р. на первом графике 
рис. 2 равно значению на графике, умноженному на коэффициент к. 

Ро - 


1-й способ 2-й способ 


р. к=1 р. к=1 
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Рис. 2. Верификация ГУ для цилиндра при 8-ми элементах базиса 
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Эта тенденция сохраняется и при увеличении числа используемых элементов базиса. Для 61-го элемен- 


та верификация граничных условий представлена на рис. 3 (приведено сравнение силы р. на участке границы 
5. ). Если оценивать погрешность как максимальное отклонение полученной величины от заданной, то во вто- 


ром способе погрешность меньше. 


1-й способ 2-й способ 


р, к=10” р., к=10” 















































Рис. 3. Верификация ГУ для цилиндра при 61-м элементе базиса 


Далее исследуем точность решения для трансверсально-изотропного тела вращения неканонической 
формы (рис. 4). Граничные условия: 
Р=О, 5, ^5,; 
р. =0,р. =0,25, 5, | р ыы ВЫ 


р. =0,р, =-1№ 5. 2=-10<и<0,5. 





Рис. 4. Граничные условия для тела вращения 


Удержано 15 элементов базиса. На рис 5. представлено сравнение граничных условий для каждого спо- 


соба (показаны не все участки границы и компоненты вектора сил). 
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Рис. 5. Верификация ГУ для тела вращения 


Как видно из графиков, для тела более сложной формы отличие в скорости сходимости наблюдается в 
пользу второго способа. Рассмотрим задачи с участием массовых сил. Последовательность формирования ре- 
шения следующая: 

|. Задается зависимость вектора перемещения плоского вспомогательного состояния от координат 


уз р и на его основе определяется вектор перемещения пространственного осесимметричного состояния. 


2. Для такого вектора определяются параметры: 

® тензор деформаций по соотношению Коши; 

е тензор напряжений из закона Гука; 

® усилия на поверхности тела из фундаментального соотношения Коши; 

е массовые силы из уравнения равновесия. 

3. Строится точное частное решение задачи, соответствующее заданной в каждой точке тела функции 


перемещения. 


Механика 


21 


БИр://уезыиК.Аоп$а.га 


ра 


Вестник Донского государственного технического университета. 2020. Т. 20, № 1. С. 15-24. 155М 1992-5950 е155М 1992-6006 
Тези! о{ оп ие Тесйтса От»егхйу. 2020. Иа. 20, по. 1, рр. 15—24. 155М 1992-5980 е155М 1992-6006 





4. Перебирая &+В<и (и=1,2,3...), строится множество точных частных решений задачи линейной 


теории упругости для параметров: 


»® векторы перемещения и, ; 

® тензоры деформаций =, ; 

® тензоры напряжений Г, ; 

® векторы поверхностных сил р, =п-Г,; 
» векторы массовых сил Х,. 


5. Формируются  базисы пространств внутренних Во в и граничных 


Ро = о состояний, в которых соблюдаются равенства: 
6, ={,,&,,Г}, 


ИИ У 
Ух — и рьХь} . 
6. Оставляем среди этих решений только линейно независимые и осуществляем их ортогонализацию в 


соответствии со скалярными произведениями в базисах внутренних и граничных состояний: 
Г @ © 
(Вы) =[=; о; АГ 5 
7 
= (Л (1). () 
(уу, )= [рФиФа$+ | ХФиФау 
5 | 


(индексы Ги /, отвечающие за номера элементов, помещены наверх и заключены в скобки). 
7. В результате получаем базис, по которому соответствующие векторы или тензоры разлагаются в 
ряды (3) с одинаковыми коэффициентами: 


с, о 


где АХ — вектор заданных массовых сил. 
Исследуем возможность построения упругого поля при наличии массовых сил, используя соотношения 
(6) и (5). Рассмотрим первую основную задачу с неуравновешенными усилиями для трансверсально- 
изотропного цилиндра (рис. 6). Граничные условия: 
Р=О, 5 95. ; 
р. =0,р. =", 5,| 2=2, О<г<1. 


Ре. 





Рис. 6. Граничные условия для цилиндра 


Выражения для ортонормированного базисного набора компонент вектора перемещения и = {и, и} 


приведены в табл. 1. 
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Таблица 1 
Ортонормированный базисный набор компонент вектора перемещения 


ПО О ОИ У 
О О ОО 
м 


и 
) 


2 
0,0908 -0,0454? —0,4347=? —0,0511=' 
—0,17577= +0, 04392? —0, 01957? +0,3661=? —0,04572* 

—Н 0,288172 —0,0804727 0,1447? —0,720927 +0,12067°27 +0,08372° 


к — 





Для получения строгого решения потребовалось 11 коэффициентов Фурье, ненулевые значения: 
с = 0,2178; с, =-0,1226; с, =0,2377; с, =-0,0732; с. =0,0247; с, =0,1443; с, =0,1443. Приведем выра- 
жения для перемещений и массовых сил (удержано 4 знака после запятой): 
и=0,2592"+0,0367 72—0,0696у2? —0,011672°; 


и’= 0, 01837? —1,14972+0,06967?2—0,0942? +; 
+0,01747?2? +0,09672° +0,01212“; 
В =-0,2814"-—0,140772; Я =1,2012-0,257* —3,60382—0,92°. 
Если построить базис внутренних состояний через плоские вспомогательные состояния, образованные с 
помощью монома 2”у’, то ортонормированный базис и коэффициенты Фурье изменятся. В этом случае реше- 


ние окажется тоже строгим и примет вид: 
и= 0, 0363+—0,091372+0,03097°2—0,0287 72°; 
и’= 0,04567* —0,00777* —0,41482- 0,0287 7*2+0,025727 +0,03992°; 
В =-2,0145"-—1,617872; Й =-0,257* —1,48932. 
Аналогично можно получить другие частные решения задачи, используя при формировании базиса по- 
линомы разного вида, например, 2’у’ +2” и др. 


Исследуемый подход позволяет получить множество решений одной краевой задачи теории упругости 
при наличии массовых сил. Подвержены перемещения и массовые силы, сочетания которых дают распределе- 
ние напряжений, удовлетворяющих на границе заданным силам. Первый способ лишен этой особенности, т. к. в 
нем массовые силы входят в состав заданных условий и задача состоит только в отыскании поля перемещений. 

Обсуждение и заключения. Второй способ решения задачи обладает наилучшей сходимостью. Кроме 
того, в отличие от первого способа при вычислении скалярных произведений в процессе ортогонализации и при 
определении коэффициентов Фурье, во втором способе не используются деформации и перемещения. Здесь 
формируется базисный набор напряжений и его след на границе — базисный набор поверхностных сил, с по- 
мощью которых проводится ортогонализация и построение рядов. Это означает, что при вычислении скаляр- 
ных произведений не возникает погрешность, связанная с составляющими, отвечающими за жесткое переме- 
щение, которое может возникать при формировании базиса [11]. 

Если рассматриваются задачи с массовыми силами, то второй способ может оказаться полезным при 
формировании множества частных решений, напряжения которых удовлетворяют определенным условиям на 
границе. Эти решения могут быть использованы в качестве базисных в более сложной задаче, а также быть по- 
лезными при определении упругих полей, реализуемых от фиктивных нагрузок, возникающих в результате 
применения метода Пуанкаре [7, 8]. 
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Проанализирована точность решения задач теории упругости методом граничных состояний при ис- 
пользовании разных подходов в построении скалярных произведений. Решение задачи линейной теории упру- 
гости, использующее представление общего решения уравнения Ламе в виде ряда Фурье по базисным функци- 
ям и выражение (6) в качестве ортогонализатора этих функций, имело наилучшую сходимость. 
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